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Pulse splitting is a crucial and common process in nonlinear fiber optics. When an intense 
laser pulse is launched into a highly nonlinear fiber, a stream of fundamental solitons is generated, 
their temporal separations increasing during propagation. This is due to the onset of a variety 
of perturbations, including higher-order dispersion and the Raman effect. Recently, it has been 
experimentally observed that the well-known law determining the amplitudes and the temporal 
widths of each soliton, however, breaks down due to the unexpected formation of metastable 2- 
peak localised states with constant temporal separation between the two maxima. In the vicinity of 
certain 'magic' input powers the formation of 2-peak states is quite common in many types of highly 
nonlinear photonic crystal fibers. In this study, we provide a full theoretical understanding of the 
above recent observations. Based on a 'gravity- like' potential approach we derive simple equations 
^^' for the 'magic' peak power ratio and the temporal separation between pulses forming these 2-peak 

^ \ states. We develop a model to calculate the magic input power of the input pulse around which the 

phenomenon can be observed. We also predict the existence of exotic multipeak states that strongly 
violate the perturbative pulse splitting law, and we study their stability and excitation conditions. 
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Pulse splitting is a well-known process occurring in the very initial moments of pulse propagation in nonlinear 
optical fibers [T, z]. According to the most accredited theory of pulse splitting [1, 2], in the femtosecond regime, 

O \ higher-order solitons are affected by stimulated Raman scattering (SRS) and higher-order dispersion terms, becoming 
unstable and eventually breaking up into several fundamental solitons. Explicit expressions for the peak power Pj 

Qh| and temporal width Tj of the j-th fundamental soliton created in the splitting process are given by (1 < j < N): 

V5 ■ Pi= Po(2N - 2j + l)ViV^ (1) 

.O; Tj=To/{2N-2j + l), (2) 

^,^. where Pq and Tq are the peak power and temporal width of the initial hyperbolic secant pulse and TV is the so-called 

r"| ■ soliton order, see Refs. [l|,l2|- After splitting, intrapulse SRS causes fundamental solitons to shift continuously to lower 

O i' frequencies via the Raman soliton self- frequency shift (RSFS) [3|, l4i]- The rate of RSFS is proportional to 1/T^ 0, S]- 

Thus, pulse splitting together with RSFS eventually leads to the complete breakup of a higher-order soliton, which 

ejects a stream of fundamental solitons one after another. According to Eqs. ([T][2]), solitons that are ejected earlier 

^ . have higher amplitudes, shorter duration, and, as a result, demonstrate stronger RSFS. This leads to a sequence of 

C^ ' fundamental solitons with different carrier frequencies, which constantly increase their temporal separations during 

CN . propagation due to the RSFS-induced walk-off effect. 

Q^ ' However, this well-known description ignores some processes that could affect the soliton dynamics, such as inter- 

^^ , action between solitons, and the emission of dispersive wave radiation near a zero group-velocity dispersion (GVD) 
l/^ ■ point. In a series of recent experiments (see Refs. |5, 6]), short and intense pulses have been launched in highly non- 
^^ \ linear photonic crystal fibers (PCFs, see also Ref. [7.]). The unexpected formation of long-lived, 2-peak soliton states 
^^ ■ for specific 'magic' input pulse energies was observed [5|, [6[. In these experiments, it was clear that this phenomenon 
is quite universal for a large variety of highly nonlinear optical fibers, especially in the deep anomalous GVD regime 
far from zero GVD points. Such localised states always show peaks of differing amplitude with a well-defined peak 
power ratio (r c^ 0.73 in the mentioned experiments). Moreover, the temporal separation between the two peaks is 
^ . also uniquely determined by r and the Raman parameter tr. Thus, under certain circumstances a pair of solitons 
^ ■ with unequal amplitudes can form a localised state, which can propagate over long distances even in the presence of 
Raman effect. 

Such soliton states were numerically discovered and studied in 1996 by Akhmediev et a/. [8]. Single peak and 2-peak 
states were found by solving the simplified nonlinear Schrodinger equation (NLSE), see Eqs. (O - dU below. The 
analysis in Ref. [8] is limited to two aspects: (i) finding the localized 2-peak solutions of Eq. (|4|), and (ii) the question 
of numerical stability of such solutions. 

In this paper, inspired by the recent experimental observation of the Raman multipeak states in Refs. [5|, Q and 
by the early numerical findings in Ref. [8], we investigate the issue of Raman multipeak states in detail. We propose 
qualitative and quantitative explanations for all the experimental observations, based on the powerful concept of 
gravity-like potential introduced in Refs. |8l-[lH. The general point-of-view introduced by our approach allows one 
to make several important predictions to be tested in future experiments, namely the possibility of forming Raman 
soliton states with more than two peaks in PCFs, thus leading to a complete violation of the soliton splitting law of 
Eqs. ([T]|2|), and to new ways to manipulate SCG in microstructured fibers by controlling exotic states of light in the 
fiber. 



RESULTS 



A. Multipeak solitons 



A simplified form of NLSE can be written, in dimensionless units, as 
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where t/j is the electric field envelope rescaled with the fundamental soliton power Ps = 1/^2 1/(7^0)5 /^2 being the GVD 
coefficient at the reference frequency. Variables z and t are dimensionless space and time, respectively rescaled with 
the second order dispersion length Ld2 = ^o/l/^2| and with the input pulse duration to. The last term of Eq. (jSj, 
responsible for the Raman effect, produces a constant RSFS for solitons in silica fibers [3|, 13, tr = Tr/Iq being a 
small parameter, where Tr is the Raman response time, approximately equal to 2 fs in silica. Solitons subject to the 
Raman effect shift continuously towards the red part of the spectrum due to RSFS, leading to a constant acceleration 
of the soliton in the time domain [3, 4]. In the reference frame of an intense accelerating soliton, one can operate the 
Gagnon-Belanger transformation '0(^,t) = /(^)exp \iz [q — 6^z^/3 + 6t)], where ^ = t — hz'^ /2 and h = ?>2TRq^ /lb. 
This leads to the following (in general complex) ODE for /(^) |8|, i9|]: 



'-Mi-{q + bOf+\f?f-rRf{\fW = Q, 



(4) 



where q is the wavenumber of the strongest soliton, proportional to its peak power. If one neglects all the nonlinear 
terms, Eq.(|4|) would correspond to the stationary Schrodinger equation for a unitary mass particle of energy q subject 
to a 'gravitational' potential U{^) = b^ [IQ, lU. Such gravity-like behavior of the Raman effect in fibers has recently 
emerged as one of the most important effects in the development of the high-frequency part of the supercontinuum 
generation (SCG), see Refs. IM^- 
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FIG. 1: Profiles of multipeak soliton solutions. Blue dashed lines show the 'gravity- like' potential U{^) created by the 
most intense soliton in its leading edge into which all other solitons with decreasing powers can be progressively fitted in an 
'organ-pipe' fashion. 

Numerical solitary- wave solutions of Eq. (j4]) with one and two peaks were found by Akhmediev for real / [8] . For 
any value of the soliton wavenumber q, such solutions always show an Airy tail on the leading edge of the pulse, 
the peak power ratio r and temporal separation ^0 of the maxima in 2-peak solutions being mainly set by material 
parameters, in particular by tr. In Ref. [8|, moreover, it was established that these Raman localised states are 
metastable, propagating for many tens of dispersion lengths before eventually disappearing due to internal collapse. 

In our investigation, we have found a variety of real solutions of Eq. (J4]) with an arbitrary number of peaks. In Fig. 
[T]we show the temporal profiles of 1-, 2-, 3- and 4-peak solitons found by solving the boundary value problem (BVP) 
by means of a shooting method with appropriate boundary conditions, for tr = 0.1. Even though these solutions have 
Airy tails on the leading edge of pulses (due to tunneling of the solutions of the linearized Schrodinger equation), such 
tails are quite small when using the physically relevant parameters, see Fig. [H Airy tails will be more pronounced if 
the slope b (or parameter q which is proportional to \/b) of the gravity-like potential gets larger. The construction 



of multipeak solutions of Eq. (j4]) can be understood by looking at the potential term q -\- b^ exhibited by Eq. (j4]). 
In fact, once the peak power {2q) and the position ^ = ^i of the strongest soliton are given, the second soliton, with 
peak power 2rq (r G [0, 1]) and position ^ = ^2, can be arranged only at a temporal separation [see Fig. [Hb)]: 



Co = a - 6 ^ q{l - r)/b. 



(5) 



This simple approximate geometrical reasoning allows one to find a relation between the temporal separation ^0 and 
the magic peak power ratio r. It is important to note that one can in principle arrange any number of individual 
solitons into a multipeak state, in an 'organ-pipe' fashion, if the intensity of the strongest soliton is large enough. 
The strongest soliton creates the gravity-like potential in its leading edge into which all the subsequent solitons, with 
progressively decreasing peak powers, can be fitted according to the qualitative geometrical relation of Eq. (J5j), see 

Figm 
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FIG. 2: Temporal evolution of multipeak solitons. 

The temporal evolution of 1-, 2-, 3-, and 4- peak solitons is shown in Fig. [2] (a,b,c,d), respectively, where the 
corresponding profiles shown in Fig. [T]are used as the initial input values for simulation of Eq.(j3]). Apart from the 
single peak soliton, all other multipeak solutions are not stable, and after some propagation length they will collapse, 
as clearly shown in Fig. [2j The more peaks these localised states have, the less robust they are, and the smaller the 
propagation length required for them to collapse. In contrast to multipeak solitons, the 1-peak soliton is quite stable, 
its profile (including its amplitude and width) being stationary to a high degree of approximation, as shown in Fig. 
Ha). 

B. Criteria of 2-peak state formation: magic peak power ratio and temporal separation 



In order to understand the underlying physics that allows the formation of the above Raman localised states we 
analyze the 2-peak soliton using a Lagrangian method. Two mechanisms are involved in the evolution of 2-peak 
solitons: (i) the intersoliton forces (ISFs) of two single initially in-phase solitons in the absence of the Raman effect 
[l2|,[l3|; and (ii) the increasing separation of these two single unequal- amplitude solitons due to RSFS in time domain. 
Here we treat these two mechanisms independently and require that they have to compensate each other in order to 
obtain a localised state. Following a theory formulated in Refs. [Iz, 13], due to ISFs two initially in-phase solitons of 
amplitudes ^/2q and y^2qr with initial temporal separation ^0 in the absence of the Raman effect will pulsate with an 
oscillation period equal to Lpi 



Ln 



2{2uy 



■ exp 



(6) 



where u = y/2q{l + v^)/4, yo = 4i^^o- The two solitons have been assumed to have a vanishing phase difference (apart 
from a constant global phase across the pulse, that can be assumed to be zero due to the gauge invariance of the 
complex NLSE), since exact numerical soliton solutions of Eq. (H]) are always purely real. On the other hand, due to 
the Raman effect, each single soliton after a propagation length Lp will acquire some additional temporal separation 
ti^2 with respect to its own initial position: ti = \bL'^ = ^TRq^L'^] ^2 = ^TRr'^q^L'^. Thus, in order to obtain the 
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FIG. 3: (a) Plot of the peak poAver ratio r and (b) the initial temporal separation ^o as a function of q. Blue dots 
and red solid lines indicate the numerical and analytical results, respectively. Parameter tr = 0.1. 

soliton state where the temporal separation of two single solitons is maintained during propagation, the following 
condition must hold true: ti — ^2 = ^o, or: 
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This condition together with Eq. (J5]) leads to a transcendental equation for the peak power ratio r: 
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From Eqs. ([5j) and (J8]) one can easily calculate the peak power ratio r and the temporal separation ^o for a given 
wavenumber q of the strongest soliton. Figure [3] shows the dependence of the peak power ratio r and the temporal 
separation ^o on the parameter q, where red sohd lines indicate the analytical (approximate) results, obtained from 
From Eqs. (|5]) and (|8]), whereas the blue dots show the numerical (exact) results, obtained by the shooting method 
for solving Eq. (|4]). It is clear from Fig. [3] that the approximate and exact results are in good qualitative agreement, 
especially when the wavenumber q of the strongest soliton gets smaller. This behavior is expected, because in this 
case the slope of the potential U{^) shown for the parameter b ^ q^ will also get smaller, so that the Airy tail will be 
negligible. Under these circumstances, the ensemble of two solitons that we use to analytically obtain Eqs. (J6j) - (J8j) 
will be closer to the 2-peak numerical solutions with smaller Airy tails. 



C. Excitation of 2-peak states 



In this Subsection we study the excitation conditions for Raman localised states by numerically integrating Eq. 
(|3j) . As input condition for Eq. (J3j) we use the combination of two hyperbolic secant pulses with an initial frequency 
detuning a as follows: ^input(^) = \^sech{^/2qt) ex.p{—iat) + ^/2qrsech[^/2qr{t — ^o)]- We launch the pulses for a 
dimensionless propagation length z = 160. For each value of the peak power ratio r we need to find the corresponding 
initial temporal separation ^o such that the temporal separation of the two pulses is as stable as possible during 
propagation. The results of this numerical analysis are depicted by the green curves with square markers in Fig. 
|4l^a,b,c,d) for four values of q = 0.369, 0.3164, 0.1151 and 0.0616, respectively, for the case when the two solitons have 
initially the same frequency (a = 0). The big black dots (we refer to these special points (r, ^o) as the magic points) 
shown in Fig. |4] represent the exact numerical soliton solutions of Eq. (j4]), see also Fig. [IJb). 

The red straight lines in Fig. |4]are calculated based on Eq. (|5]), and the blue solid curves are obtained from Eq. 
([7]). The crossing points of these curves correspond to the solution of r as found by our qualitative model of Eq. (J8]). 
The discrepancy between the analytical curves (blue solid ones) and the numerical curves (green curves with square 
markers) gets smaller for larger r and ^o- 

When the parameter q is large enough (see the green curves with square markers in Fig. |4l^a,b) with q = 0.369 
and 0.3164, for propagation length z = 160 and initial frequency detuning a = 0) we do not find any solution (r, 
^o) very close or below the magic points (single big black dots in Fig. |4]) so that the temporal separation of the two 
pulses is preserved during propagation. Of course, the localised states can survive for shorter propagation lengths 
(e.g. z 0:^ 10), but they collapse during further propagation. This behavior takes place even exactly at the magic 
points. This means that the Raman localised states studied here (including magic points) are not mathematically 
stable, but metastable as found in Ref. [8|. Nevertheless, these metastable states can propagate for surprisingly long 
distances, as experimentally observed in Refs. [5|, Q. If we go upwards along the green curves with square markers 
in Fig. m the localised states survive for longer propagation lengths. This is expected, because when the peak power 
ratio r and the temporal separation ^o become larger, the amplitude difference between the two pulses falls. As a 
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FIG. 4: Peak power ratio r versus initial temporal separation ^o. (a,b,c,d) Results for q — 0.369, 0.3164, 0.1151 and 
0.0616, respectively. Red dot lines are calculated based on Eq. (|5]), blue solid curves are obtained from Eq. (|7|), crossing points 
of these curves correspond to the solution of Eq. (|8]). Single big black dots are magic points. Green curves with square markers 
are obtained by numerically modeling Eq. (|3|) for a propagation length 2;= 160 with two sech solitons initially having the same 
frequency [a — 0) as initial input values. The curve with triangular markers in (b) shows the case a — —0.2. 

result, the Raman effect is almost the same for them, they are located further from each other, and thus their ISFs 
get smaller. In the extreme case when two pulses have the same amplitude and are located infinitely far from each 
other, this localised state will obviously survive for quite a long propagation length, because each fundamental soliton 
is very robust [8]. 

The curve with red triangular markers in Fig. IH^b) shows the case when the two pulses have initially different 
frequencies (the initial frequency detuning a = —0.2). In this case, for the same peak power ratio r, the initial 
temporal separation ^0 between the two pulses is smaller compared to the case when a = 0. But during propagation 
the temporal separation gets larger and asymptotically stabilizes. In the case a = the temporal separation is almost 
the same during the whole propagation length. 



D. Magic input powers 



Once the values of r and ^0 have been established, one can explicitly calculate the 'magic input power' of an intense 
input pulse around which the formation of 2-peak localised states will appear during the pulse splitting process. 
According to the accepted theory of pulse splitting under the effects of Raman perturbation [ll, |2| , an intense A^th- 
order soliton V^input = A^sech (t) will split into A^ individual fundamental solitons of form y^sech ( y^^t) , where 
the powers Pj are given in Eq. ([T]). The condition for which two successive solitons (for instance j and j + 1) have 
a peak power ratio equal to r, turns into the equation Pjj^i/Pj = r, which has the physically relevant solutions 
A^mag,j = (j — 1/2) + [1 — aA"]"^- The first two solitons with stronger peak powers (j = 1 and j — 2) will have the 
correct amplitude ratio for A^mag,i = 1/2 + [1 — \f^\~^ - For the case r = R^ c^ 0.73, one obtains Amag,i — 7.4 which 
corresponds surprisingly well to the measured value (A" c^ 8), see Refs. [5|, |6|. According to the above straightforward 
theory, only one 2-peak soliton can be created starting with one intense input pulse, since at each magic input power 
only one pair of solitons will possess the magic peak power ratio r. It is also clear from the above relation that, for 
a given value of the peak power ratio r, a higher input soliton order A^ is required to form a localised state between 
the j-th and the (j + l)-th fundamental solitons with higher value of j. 

Figure [5] shows the dependence of the spectra at the output for z = 0.45 as a function of the input soliton order 
A' when using a soliton of the form V^input = A'sech (t) as the initial value for numerical simulation of the generalized 
NLSE (GNLSE) [Eq. (|9|)] with the full convolution for the nonlinear term: 






Rit')\i>iz,t-t')\'^dt' = 0, 



(9) 



where the response function R{t) includes both the electronic (Kerr) and vibrational (Raman) contributions [1^. One 

has R{t) = (1 — 0)5{t) + Oh{t)^ where h{t) = [rf + r2][rir2]~"^ exp(— t/r2) sin(t/ri) is the Raman response function 
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FIG. 5: Spectral evolution with increasing soliton number N. The propagation length is z = 0.45. 

of silica, with ri = 12.2 fs and T2 = 32 fs; S{t) is the Dirac delta function. Coefficient parameterizes the relative 
importance between Raman (non-instantaneous) and Kerr (instantaneous) effects, and for silica :^ 0.18. Parameter 
tr in Eqs. (J3j) - (j4]) is related to ri^2 as follows: tr = 2^rf r2/[to(rf + rf )], to being the input pulse duration. With 
suitable approximations Eq. (|9]) can be shown to reduce to Eq. ([3]). In Refs. [H, [6| the numerical pulse propagation 
analysis includes higher order dispersion terms. In this paper, as seen from Eq. (|9]), we limit ourselves to second order 
dispersion because we wish to eliminate the phenomenon of resonant radiation emission from solitons [15], which was 
present in the experiments of Refs. [5|, l6|, but was not shown there. Operating far from any zero-GVD point increases 
the chance of observing Raman localised states, since their formation will not be disturbed by the small amplitude 
background waves generated in the normal dispersion regime close to a zero-GVD point. However, many solid-core 
PCFs have a zero-GVD in the near-visible and a broad region of anomalous GVD in the far infrared. In such fibers, 
even though solitons emit radiation if launched near the zero-GVD point, they shift continuously towards the deep 
anomalous GVD region, where in many cases the GVD is nearly constant for a broad range of frequencies - hence 
the constant GVD approximation used in our model of Eq. (|9]). 
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FIG. 6: Formation of 2-peak localised states. 

12.075, 12.65, and 14.65, respectively 



(a,b,c,d) Temporal evolution of pulses with input soliton order A^ = 9.63, 



Within the three regions enclosed by white ellipses in Fig. [5] soliton pairs can be formed between adjacent funda- 
mental solitons ejected during the break-up of the input soliton. Four values of soliton order N = 9.63, 12.075, 12.65, 
and 14.65 are used in these regions to show the pulse temporal evolution in Fig. [6fa,b,c,d), respectively. The localised 
pair formed by the 1st and the 2nd fundamental solitons is depicted in Fig. [6l[a) [N = 9.63 is located in the left-hand 
ellipse in Fig. (J5j)], the inset showing the details of this localised pair and the third separate fundamental soliton. 
In this case, after a propagation length z = 1, the peak power ratio of two pulses forming the Raman state is r = 
0.7534, corresponding to the calculated value A/'^ag = 8.07. Here we should mention that during pulse propagation 



the two pulses forming the pair wih exchange energy, so that the peak power ratio r wih also slightly change during 
propagation. Even though the model used in this section to calculate the magic power is quite crude and qualitative, 
the correspondence with the experiments is surprisingly good. 

Figure [6l^b,c) with input soliton order N — 12.075 and 12.65 (the middle ellipse in Fig. [5]) shows the states formed 
by 2nd with 3rd, and 3rd with 4th fundamental solitons, respectively. The pair formed by 4th and 5th fundamental 
solitons is depicted in Fig. [6l^d) with input soliton order N = 14.65 (the right-hand ellipse in Fig. [5j). All the localised 
states shown in Fig. [6] are metastable, collapsing after sufficiently long propagation lengths. The results shown in Fig. 
[6] are remarkable in the sense that they demonstrate that one can generate pairs for any selected adjacent solitons, 
simply by properly adjusting the input pulse energy. 

In case of solitons represented by magic points, see Figs. [1] and [21 not only the temporal separation ^q^ but also 
the intensity profiles of each pulse, and thus, the peak power ratio r, will be preserved where they are stable. So, 
the localised states shown in Fig. [6] are not solitons presented by magic points. Moreover, solitons corresponding to 
magic points are solutions of the simplified model of Eq. (J3]), not of the more complete equation (|9]). 
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FIG. 7: Formation of a 3-peak Raman state. 

Intensity profile of this 3-peak Raman state at z - 



(a) Temporal evolution of a pulse with input soliton order N = 13.81. (b) 
0.8. 



So far we have only investigated the 2-peak states of the GNLSE [Eq. ([9|)]. Apart from 2-peak Raman states, the 
simplified form of NLSE [Eqs. (J3]|4])] also possesses localized solutions with an arbitrary number of peaks [provided 
that there is enough 'space' to arrange them according to Eq. (|5])], as shown in Fig. [H It turns out that the GNLSE 
[Eq. ([9])] also shows evidence of formation of 3-peak states in its dynamics, contrary to the oversimplified conclusions 
that can be made by using Eqs. ([l]|2]). This leads to the important prediction that in properly designed realistic PCFs 
there must be a strong violation of the soliton splitting law given by Eqs. (PE]). Indeed, Fig. [Tl^a) shows the temporal 
evolution of an intense pulse with the input soliton order N = 13.81 (the right-hand ellipse in Fig. [5j). During the 
pulse splitting process, a 3-peak state is formed by the 2nd, 3rd and 4th fundamental solitons. The intensity profile 
of this 3-peak state at z = 0.8 is shown in Fig. [Tl^b). Again, this intensity profile demonstrates that solitons organize 
themselves in an organ-pipe fashion, as consistent with the gravity-like potential approach used here, see also Fig. 
He). 
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FIG. 8: Experimental data. Power dependence of the spectra obtained experimentally at the end of a highly nonlinear PCF 
(18 cm long). 



Two-peak Raman states were experimentally observed and reported in Refs. [5|, [6[. From the same experimental 
data we have also observed some indication of the existence of 3-peak states. Figure [S] shows the power dependence 



of the spectra obtained experimentally at the end of a highly nonlinear PCF (18 cm long) [sl, [6|. At ~ 100 mW, a 
2-peak Raman state is formed by the 2nd and the 3rd fundamental solitons, where at ~ 200 mW one can see some 
signature of a 3-peak state formed by the 1st, 2nd and 3rd fundamental solitons. These three solitons are close to 
each other in frequency. At the same time, the presence of strong spectral fringes at ~ 200 mW (not shown here) 
shows that they are also close in time. As mentioned above, 2-peak states are more robust than 3-peak states, but 
the existence of solitonic states with more than 3 peaks would unambiguously confirm our theoretical model. 

II. DISCUSSION AND CONCLUSIONS 

The formation of multipeak localised structures supported by the Raman effect turns out to be a fundamental 
and surprisingly common process occurring in the region of deep anomalous dispersion in photonic crystal fibers. 
The 'gravity- like' potential approach turns out to be a fruitful tool, leading to a simple relationship between two 
important parameters of 2-peak states, namely the peak power ratio r and the temporal separation ^o between 
two pulses forming the pair. Our approach explains in a natural way the 'organ-pipe' arrangement of multipeak 
states, where the strongest soliton creates a potential in its leading edge and subsequent solitons with monotonously 
decreasing peak powers can be fitted within this potential, see Fig. [H The formation of pairs can be explained 
through the following simple mechanism: the RSFS induces an increasing separation between two unequal-amplitude 
solitons, that can only be counterbalanced by intersoliton forces. Analyzing these two effects together with the results 
obtained from the potential approach led us to a simple transcendental equation for the peak power ratio r, and thus 
also for the temporal separation ^o- This transcendental equation gives results that are in good agreement with the 
exact numerical solutions. The physical interpretation of 2-peak state formation in our model is thus quite simple 
and intuitive. 

We have used a qualitative model to predict the magic pulse power needed to generate a 2-peak state during the 
pulse splitting process. The predicted values are again in good agreement with experimentally measured ones. In 
addition to 2-peak states, we have also investigated n-peak states with n > 2. We have demonstrated numerically - 
using some hints from experiments - that the complexity of the GNLSE allows the existence of 3-peak state. Even 
though such states are much less robust than 2-peak states, they can survive for at least several tens of dispersion 
lengths. It is hoped that such 3-peak states can also be experimentally observed in an unambiguous way in the near 
future. 

We expect that the results reported here, apart from their obvious fundamental interest, will allow researchers to 
manipulate more effectively the dynamics of solitary waves and supercontinuum generation in solid-core microstruc- 
tured fibers with a properly designed group velocity dispersion, and in hollow-core fibers filled with Raman-active 
gases. 
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